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NATIONAL ADVISORY COMMITTEE FOE AERONAUTICS 
RESEARCH MEMORANDUM 

PRELIMINARY ANALYSIS OF PROBLEM OF DETERMINING EXPERIMENTAL 
PERFORMANCE OF AIR-COOLED TURBINE 
I - METHODS FOR DETERMINING HEAT-TRANSFER CHARACTERISTICS 
By Herman H. Ellerbrock, Jr., and Robert R. Ziemer 


SUMMARY 

Little material lias been published on the methods of investi- 
gating air-cooled turbines and analyzing, the data obtained to eval- 
uate the turbine performance, which includes the heat-transfer and 
cooling-air-flow characteristics in addition to the performance 
characteristics normally obtained in uncooled turbines. Some anal- 
ysis and study has been made by the NACA Lewis laboratory in con- 
nection with investigations that are to be made on air-cooled 
turbines. Methods that are required to determine the character- 
istics mentioned are presented. The methods have a form in which 
dependent parameters are expressed as functions of independent 
parameters, the functions generally being unknown. Methods of 
experimenting to determine these functions are suggested. Although 
the forms of the formulas are, for the most part, unchecked experi- 
mentally, it is believed that the information will be useful in the 
turbine- cooling field. 

The methods that are expected to permit determination of the 
heat-transfer characteristics of an air-cooled turbine from specific 
investigations that must be conducted are presented herein. In some 
cases 'discussions that lead to the suggested form of the formulas 
are given in order that readers unfamiliar with the subject may 
obtain a knowledge of the general background of heat transfer. 


. INTRODUCTION • 

Little material is available at present on the methods of 
investigating air-cooled turbines and analyzing the data obtained 
to evaluate their performance. The performance of turbines formerly 
included such factors as efficiency and power but did not include 
characteristics such as heat transfer and cooling-air flow because 
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the turbines were uncooled. Some study and analysis has been made 
at the NACA Lewis laboratory on the problem of investigating and 
evaluating the performance of air-cooled turbines. Because manu- 
facturers are applying cooling to experimental turbines, the results 
of the studies at the Lewis laboratory, although lacking experimental 
verification, are presented in a series of reports. 

The reports will present suggested formulas that are required 
to determine air-cooled turbine performance. The formulas will be 
set up so that certain dependent parameters are functions of inde- 
pendent parameters, the functions in general being unknown. Methods 
of investigation to determine these functions will be suggested. 

The suggested methods for determining the heat- transfer character- 
istics, the cooling-air flow characteristics, and the power and the 
efficiency of air-cooled turbines will be presented. These reports 
will be concerned only with the turbine. The problem of power loss 
in the engine due to bleedoff of cooling air from the compressor 
when high pressure is required and other engine problems are not 
included in the evaluation of turbine performance. 

Methods based upon present knowledge that will permit deter- 
mination of the heat-transfer characteristics of a turbine from 
specific investigations, which must be conducted, are suggested 
herein. In some cases discussions are given In order that those 
unfamiliar with the subject may obtain some knowledge of the general 
background of heat transfer. From the heat- transfer characteristics, 
factors such as the maximum permissible gas temperature, the amount 
of heat added to the coolant, and the cooling-air requirements, which 
are of primary importance in determining the power of the turbine or 
the thrust of the engine, can be obtained for conditions other than 
those investigated. 


METHODS OF EVALUATING HEAT-TRANSFER CHARACTERISTICS 
Factors Required 

It is important to be able to determine the gas temperature at 
which a cooled turbine operates safely for a given coolant flow and 
coolant temperature (or, the coolant conditions that are required 
to operate at a given gas temperature without exceeding blade-stress 
limitations). From such a determination, the net performance of the 
turbine or engine can be found by use of suitable formulas. 

In order to determine the safe operating conditions, the 
allowable-blade- temperature- distribution curve and the actual-blade- 
temper ature- distribution curve must be found. The actual-blade- 
temperature- distribution equation involves heat-transfer coefficients 
and effective fluid temperatures for which formulas must be available. 
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As heat is absorbed by the coolant, the performance of the 
turbine varies. Consequently, a parameter involving the heat to the 
coolant enters into the performance determination and formulas for 
the calculation of this parameter must be set up using test data. 
Discussions of the details of the heat-transfer formulas and sec- 
ondary formulas, which are required to obtain the heat- transfer 
formulas, are given. The discussions are concerned with the rotor, 
except where otherwise noted. The last part of this section is 
concerned with the cooled nozzle for which the methods of evaluating 
the rotor heat- transfer characteristics can be applied with only 
small changes. 


Allowable Badlal Blade Temperature Distribution 

The manner of determining the coolant flow required for safe 
operation of the turbine at given effective gas and coolant tem- 
peratures is shown in figure 1. For each turbine speed an allowable- 
radial- temperature-distribution curve exists. Also for a given 
effective gas temperature Tg^ e and effective coolant temperature 

T a ^ e , the coolant flow determines the actual radial temperature 

distribution, as illustrated in figure 1. A particular coolant flow 
gives a curve such as A (fig. l), which will contact the allowable 
temperature curve at some point such as C. This point then rep- 
resents the lowest allowable coolant flow for the given operating 
conditions. A lower flow would cause curve A to move upward and 
exceed the allowable temperature over a portion of the blade, causing 
blade failure. 

The allowable blade temperature depends on the stress character- 
istics of the blade. A brief discussion of this subject is given in 
reference 1. The conclusion of the discussion is that, for the 
present, the allowable blade temperature distribution must be deter- 
mined from the radial centrifugal stress distribution in the blade 
and from stress-rupture data. The method of determining this allow- 
able blade temperature distribution, which is quite simple, is 
reviewed herein for convenience. 

The centrifugal-stress distribution is determined by calculating 
the stress, increment by increment, from the tip to the root of 
blade. (See fig. 2(a).) The centrifugal force of an increment is 

centrifugal force = A B ^ aT dr CO r (1) 
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The area A B,av is the "blade cross-sectional area at the center of 
the increment in question. (All symbols are defined in appendix A. ) 
The blade stress at any section is the sum of the centrifugal forces 
acting on all the increments from the section to the tip divided by 
the area of the section, or 


S 


B + p_ A dr (0 2 r 
B B,av 

144 A b 


( 2 ) 


vhere A B ^ av , A B , r, and dr correspond to the increment under 

consideration and B is the sum of the centrifugal forces on pre- 
ceding increments. The area A B is that of the base of the incre- 
ment. (See fig. 2(b).) The stress at any portion of the blade can 
be determined by this step-by-step method. In some cases where the 
areas of the blade sections can be represented as a function of r, 
the stress distribution can be obtained by integration. 


"When the centrifugal- stress distribution is known, the allowable 
radial temperature distribution can be determined from stress-rupture 
data. The 10,000-hour life stress-rupture curve has usually been used 
in order to obtain a factor of safety because of the vibratory and 
thermal stresses that are neglected. 


Actual Blade Temperature Distribution 

The actual blade temperature distribution formula can be con- 
siderably simplified if the following assumptions are made: 

(a) Eadiation effects are negligible. 

(b) Wall thickness is negligible insofar as heat transfer is 
concerned ; that is, no temperature drop occurs through the wall. 

(c) Conduction of heat in the radial and peripheral directions 
in the blade wall is negligible. 

(d) Blade mean cross section applies to the entire blade when 
average-value calculations are made. 

(e) Local convection heat-transfer coefficients on the Inside 
and outside of the blade are constant spanwise and peripherally at 
given operating conditions. 
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(f) Local effective combustion- gas temperature (temperature 
effecting beat transfer on outside of blade) is constant over 
outside blade surface. 


With these assumptions the equation for the blade temperature at 
any point x along the span from heat-transfer considerations is 


T «.e - Tb.i 

T -T 
g,e a,e,h 



2 . 
0) v a b 


jgH i T¥ r* 

o o g,e a,e,h J 


I “ X+Y 

b 


(3) 

Equation (3) can be derived from equations given in appendix B of 
reference 2. 

Badiation between the blades and other surfaces and radial heat 
conduction has little effect on the temperature distribution except 
near the blade root (reference 2). This fact, however, does not 
appreciably affect the determination of the required coolant flow or 
the allowable gas temperature. Gaseous radiation may possibly become 
a factor to consider in some special cases such as very high gas 
temperatures and gases containing certain constituents. In general, 
however, gaseous radiation can be neglected. Assumption (b) should 
apply well for the thin- walled blades used in air-cooled turbines. 

The neglect of peripheral heat conduction, assumption (c), is war- 
ranted for the greatest part of the blade on the basis of unpublished 
investigations already made on cascades. The error involved through 
use of assumptions (d), (e), and (f) must await comparison of tem- 
peratures calculated from equation (3) and experimental values. On 
the basis of some study, it is believed that the use of these assump- 
tions will result in small error in the blade temperature distribution. 

The use of assumption (b) and neglect of heat conduction in a 
peripheral direction reduces the heat transfer to a one- dimensional 
problem, that of radial heat transfer, and equation (3) is conse- 
quently the equation for the radial temperature distribution. For 
a given set of combustion- gas and cooling-air conditions, all terms 
in equation (3) are constants except T-g x and x. Thus the 
equation reduces to one of simple form for use in the analysis of a 
turbine. In some cases the X and T terms can also be neglected 
resulting in further simplification. 

Before applying equation (3) to determine the blade temperatures 
of a given turbine, equations for BLj_, T g ^ e , and T a ^ e must 

be determined from data obtained in experiments on the turbine. 

These equations are now discussed. 
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Outside Convection Heat-Transfer Coefficient 

Method of calculation. - The heat addition to the cooling air 
as it passes through the blades, the average effective gas tem- 
peratures, and the average blade temperature must first be deter- 
mined from experimental data. With these data and the blade dimen- 
sions, the outside convection heat-transfer coefficient H Q is 
calculated. The following equation is used: 


20 = (*g, a ^ T - *B,av) (4) 

This equation again assumes that all radiation is negligible. When 
radiation is not negligible, the discussion in the previous section 
is applicable. Discussions of methods for determining T g^ e ,av>» 

and Tg^ av from experimental measurements are discussed in the 
following paragraphs. 


Heat gained by cooling air Qa. - From the general energy equa- 

tion, it can be shown that the heat gained by the cooling air in the 
blade cooling-air passage for a short distance dr can be expressed 
by the following equation: 


dq a - dfl 


c p,a 


dr. 


+ d 



( 5 ) 




«* 


where 

dq a differential quantity of heat gained by cooling air, 

(Btu/lb cooling air) 

dfi work added to cooling air by blade rotation, (Btu/lb cooling air) 
dT & rise in static temperature of cooling air, (°F) 

It can be shown that the right-hand side of equation (5) equals. 
Cp^ a dT M a where dT" a is the rise in total, temperature of the cooling 

air relative to the blade. The work is 

W^r dr 
Jg 


dO 


( 6 ) 
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Tims equation (5) reduces to 


- ep,«®"a- 


CD 2 r dr 

Jg 


( 7 ) 


By integrating from the roots to the tips of the hlad.es and hy 
assuming Cp a to he a constant hut determined at the mean of the 
static temperatures at the root and the tips, equation (7) becomes 

^a “ °p,a,av a,T “ T "a,I^ _ 2Jg (^ 2 “ 

The total heat to cooling air is then 

Qe - v aqs (9) 

Because of circulatory effects of the fluid in the passages, 
the angular velocity of the cooling air, which is the term that 
should be used in equation (8), is a little smaller than the angular 
velocity of the wheel CO, which is used in equation (8). Conse- 
quently, it is more accurate to write equation (8) as 

^a * °p,a,av ^"a,T " T "a,h^ “ f s jjjg (^T r h J 

where f s is the factor (equivalent to the slip factor of a cen- 
trifugal supercharger) to correct for the circulatory effects of 
the cooling air. The value of this factor can he determined hy 
motoring the turbine wheel at several speeds and, with no gas passing 
across the turbine blades, allow various cooling-air quantities to 
pass through the blades. Thus q ft will be 0 and f s can be cal- 
culated from equation (10) from the measured air temperatures and 
set up as a function of w a and CD. 

The total temperatures of the air relative to the blade at the 
root and the tip are determined by means of the usual formulas 
involving the recovery factor of the thermocouples, obtained from 
calibration tests, measured temperatures, and measured values of 
static and total pressures in the cooling-air passages. The static- 
and total-pressure measuring devices are placed as close to the 
thermocouples as possible. A more detailed discussion of these 
measurements is given later because obtaining measurements in 
rotating passages involves some special problems. 



8 


NACA EM E50A05 


It is important to be able to determine Qq for the turbine 
from the established formulas in order that the performance of an 
engine in vhich the turbine is to be placed can be calculated for 
any flight and engine conditions. Equation (10) can be used if 
the cooling-air total temperatures can be calculated from test- 
verified formulas. The methods of obtaining T" a for any condi- 
tion so that can be calculated are discussed in reference 3. 

Effective gas temperature e . - The temperature effecting 

heat transfer from the gas is equal to the temperature the blades 
would assume if they were thermally insulated (no heating or cool- 
ing) under the same fluid conditions for which the heat - transfer 
coefficient is being determined. The effective gas temperature 
thus is the adiabatic blade- surface temperature. The reason for 
the use of such a temperature to determine convection heat-transfer 
coefficients is discussed in reference 1. 


The adiabatic temperature of a body such as a turbine blade, 
like that for a thermocouple, can be determined, as brought out in 
reference 1, from the total and static temperatures of the fluid 
flowing around or through the body and from a recovery factor. 

Thus it is possible to conduct investigations in which no cooling 
air is used and to determine the average total and static tempera- 
tures about the blades as well as the blade temperatures. The 
recovery factor for the blades can then be found from the equation 


m” m 

A =* 1 - -Si® 

e t" _ - t 


g,av 


g,av 


( 11 ) 


where Tg^ Q is equal to the average of the measured blade tem- 
peratures T-g^ av under adiabatic conditions. 

A summary of the studies of many investigators on various 
bodies. Including that of reference 4 on turbine blades, to determine 
the parameters affecting the effective gas temperature is given in 
reference 1. In general a formula of the following form is 
applicable: 




(Pr g )“ 




( 12 ) 
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where f^ denotes an unknown function and n an exponent of the 
Prandtl number. For the case of laminar boundary layer on a flat 
plate the exponent n has been shown to be equal to l/2. 

Tests must therefore be conducted on the turbine over a range 
of Prandtl and Mach numbers at gas temperatures that are safe when 
no cooling air is passed through the blades to establish the 
function of the formula given by equation (12). 

It is possible, however, that adiabatic conditions are not 
obtainable at temperatures much different from, those surrounding 
the experimental setup. Also, over the probable range of tem- 
peratures covered in this type of investigation, the Prandtl 
number will not vary more than 5 percent and this number raised to 
a power less than 1.0 will vary even less.' The accuracy of the 
data may therefore not warrant an investigation varying the Prandtl 
number. The recovery factor Ag is calculated using equation (11) . 

After establishing equation (12), effective gas temperatures for 
use in equation (4) for the conditions of high gas temperature used 
when cooling air is passed through the blade, can be calculated 
using equation (11), the recovery factor obtained through use of the 
established equation (12), and measured fluid conditions. 


For the case of establishing equation (12) from low-temperature 
experiments with no cooling air and the case of calculation of Tg^ e 
for high- temperature experiments with cooling air to determine Hq, 


methods for determining the average relative total gas temperature 


mil 

g,av> 
number Pr 


the average static gas temperature T 


the Prandtl 


g' 


and the average Mach number 


g,av> 

M g? 3 from measurements 


made during the investigations must he known. Also the method of 
determining the average blade temperature T^^ aT in establishing 

equation (12) must be known. These methods are described in the 
following sections. 


It Is recommended that experiments with variable Mach number 
be conducted at several constant gas Eeynolds numbers to verify the 
neglect of the Eeynolds number in equation (12). (See reference 4.) 
The Mach number variation can be obtained by varying the fluid flow 
rate and turbine speed. Evidence exists that the formula for A~ 


should be based on the Mach number at the rotor-blade outlet 


%3 


(references 4 and 5). The reason for the variation of the Mach 
number basis of A_ is thought to be due to the boundary- layer 

O 


flow conditions that exist in each case. Until more data are avail- 
able, it is recommended that the Mach number that correlates the 
data he used. 
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Because the differences "between the fluid temperatures and 
blade temperatures when determining Ag are very small, the blade 

temperatures should be measured differentially with respect to 
some one fluid temperature if any degree of accuracy is to be 
obtained. 

Average blade temperature - It is possible to con- 

nect the various thermocouples used to measure the blade-wall tem- 
peratures if they are insulated from each other in a parallel 
circuit so that the reading of the measuring instrument used will 
be an average for the given test conditions imposed on the turbine. 
A circuit of this type, however, might cause mechanical difficul- 
ties when used on a rotating wheel and the average may not be an 
integrated average. 

It is suggested that individual thermocouples be placed in 
peripheral bands at various stations along the blade span. The 
thermocouple readings can then be integrated with respect to the 
distance around the blade for each peripheral band. Then if the 
integrated mean values of temperature for the peripheral bands are 
integrated across the blade span, the average blade temperature 
Tjg ?av is obtained. 

Total temperature of cooling air T" a . - The observed tem- 
peratures of the cooling air at the blade root and tip can be 
obtained by means of a thermocouple placed at each position in the 
blade cooling-air passage. It is assumed that any thermocouple 
placed at each station will be so located that a good average of 
the air temperature across the passage is obtained by the one 
reading. Because the thermocouples are attached to the blades, the 
total temperatures calculated from the observed readings will be 
values relative to the blade, which are the values required. 

It is quite likely that the thermocouple readings will be in 
error if improperly shielded and insulated because of radiation 
from the hot blade walls and conduction from the blade walls along 
the thermocouple support to the Junction. Much care and considera- 
tion should therefore be given to the installation of the 
thermocouples . 

As previously mentioned, the total temperatures T" a ^ and 
T" a ,<p required in equation (10) are calculated from the observed 

readings, the recovery factor of the thermocouple, and the pres- 
sures in the passage. The recovery factor is also a function of 
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the Mach number of the air so that the pressures at each station 
are also required to determine the Mach numbers. These pressures 
are discussed in the following section. 

Cooling-air pressures in blade passage. - The total pres- 
sures p ff a and static pressures of the cooling air at the 

blade root and tip can be obtained by means of pressure tubes 
attached to the blade and static-pressure taps in the blade wall. 

It is desirable that static-pressure tubes be located in the cooling- 
air passage in addition to use of taps in the wall in order to get a 
true average pressure. In most cases, it will probably be impossible 
to put more than one or two instruments in any one passage. If more 
than one instrument is required to obtain an average pressure, these 
instruments will have to be apportioned among several blades. 

The pressures in the rotating-blade coolant passages are 
transmitted to a point near the center of the turbine through tubing 
and then through a pressure pickup to a stationary manometer board. 
The centrifugal effects- in the interconnecting tubing must be taken 
into account when evaluating the pressures in the passage from the 
manometer-hoard readings. The equation for using the manometer- 
board readings to obtain the correct pressures in the coolant pas- 
sages is derived in appendix B. For pressures, either total or 
static, obtained with pressure tubes located at the blade root, the 
following equation is applicable: 


Rotating pressure « (manometer differential pressure ± 
pressure of room in which manometer is placed) 

2 ^a T a,t (13) 

e 

The temperature T a ^^. is the temperature of the fluid in the con- 
necting pressure tube, which is discussed in appendix B. For 
pressures at the blade tip, is substituted for in 

equation (13) . 

Average static gas temperature T g ^ av . - The ideal method of 

obtaining an accurate average static temperature would entail 
measurements through the gas passage between the blades and then 
integration of the resulting values with respect to the blade- 
passage length. This method, however, requires a large turbine and 
elaborate instrumentation. A more practical indication of the 
average static temperature can be obtained from the expression 
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T 


g>av 


T g,2 + T gj5 

2 


(14) 


vhere 



T g,a 

static 

temperature of 

T «,3 

static 

temperature of 


o 

gas behind nozzles , E 
gas downstream of rotor, ^ 


In some turbines it will probably be impractical, if not 
impossible, to place instruments at station 2. In this case, the 
temperature of the combustion gas must be calculated from conditions 
upstream of the nozzles or the condition at station 2 for specified 
static pressures will have to be calibrated with the rotor removed 
with reference to a measuring station upstream of the nozzles. 


In calculating the static conditions and velocity of the gas 
at station 2 from the measured values at station 1 for the case of 
uncooled nozzles, it is assumed that the total temperature T* 2 
and the total pressure P , g J 2 downstream of the nozzles are equal 
to the total temperature T, g i 033x1 the total pressure p* ^ 

upstream of the nozzles. The assumptions are thought to be valid 
inasmuch as only a small heat loss to the casing is neglected. If 
the nozzles are cooled, a correction of the total temperature at 
station 2 must be made to allow for the heat added to the coolant. 
The heat added to the cooling air flowing through the nozzles is 
determined from the formula 


Qs = w a,S c p , a , av (^ T, a,T,S “ T 'a,h,s) ( 15 ) 

The total temperatures of the air in the nozzle cooling passage and 
the cooling-air flow passing through the nozzles are determined 
from measurements. The total temperature of the gas at station 2 
Is then obtained from 


T 


g,2 


_ ml _ % 

8,1 Vp,s 


( 16 ) 


The specific heat of the gas at this point is determined from 
because Tg^ 2 is ^ eln S calculated. If after Tg^ 2 18 

calculated based on the average of a 11 *! ^g,2 is 
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much different from the °-p,g used at first, a recalculation of 
should be made using the c p^g based on the average tem- 
perature. With cooled nozzles, the assumption of p* g ©qual to 
p*g is thought to be more valid than before inasmuch as the 

cooled nozzles tend to stabilize the boundary layer and reduce the 
blading losses. 

The method of calculating static conditions behind the nozzles 
[knowing T' g ^ 2 > P*g,2> v g> A 2 (the nozzle-outlet area)Q 

is primarily a trihl-and- error solution. With two additional 
assumptions, however, the nozzle-outlet conditions can be quite 
closely approximated. The following equation is developed in 
reference 6 (equation (9)): 


1 



where C is the nozzle-discharge coefficient. With the nozzle- 
discharge coefficient assumed equal to 0.98, all the terms on the 
right-hand side of equation (17) are known from turbine measurements 
and the assumptions made, with the exception of the ratio of spe- 
cific heats 7 g . For the first trial, the ratio is based on T’ g g. 

After equation (17 ) has been solved for the pres sure- ratio term 
(the left-hand side of equation (17)), the pressure ratio is obtained 
by making use of one set of curves in figure 3. When the pressure 
ratio is known, the temperature ratio ^g,2/^*g,2 is obtained from 

the other set of curves in figure 3, which are based on the simple 
formula 





7 


g 


( 18 ) 
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If the value of T g ^ 2 is appreciably smaller than T' g 2 > a uew 
7 g based on the calculated T g ^ 2 is used and the procedure 

repeated. With the calculated value of T gj2 and - the value of 

obtained by normal means from measurements at station 3, the 
average static temperature can be found. (See equation (14).) 

The foregoing method applies to flows up to a Mach number of 
1.0 at the nozzle outlet and for no further expansion past this 
point. The case of supersonic flow caused by an expansion down- 
stream of the nozzles is not considered herein. 


Average total gas temperature T" g ^ av .. - The relative total 
gas temperature ^g^av can be evaluated "from the average static 
gas temperature T g ^ av and the relative average rotor-inlet and 
-outlet gas velocities, ^2,a,y and w 3,av* respectively. The 
expression used Is 


Wo + W-r _ w 

<£ , av 5,av 


T ’g,av “ T g,av + 


2 Jgc 


P;g 


(19) 


Velocity diagrams, such as shown in figure 4, are made for four or 
five sections along the rotor-blade span and values of the inlet 
and outlet velocities are integrated to obtain the averages 
required in equation (19). 


The velocity diagrams at any section are constructed as fol- 
lows: When the turbine-rotor speed and the radius to ihe section 

being considered are known, the tangential velocity u of the 
blade can be evaluated. (See fig. 4.) The angle dg Is known by 

assuming that the gas leaves the nozzles in the direction of the 
nozzle -outlet angle or by making surveys downstream of the nozzles 
with the rotor removed. The velocity Vg can be calculated using 

measured values of T g ^ 2 and - ^"g,3 ° r *>7 using those calculated 
by methods previously given and the formula 



T g,2 



( 20 ) 
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* The nozzle velocity coefficient was considered in the evaluation of 
Tg^g* The velocity Wg can then "be determined. (See fig. 4.) The 

velocity V 3 and direction angle ctg are determined hy survey 
measurements just downstream of the rotor. 

Prandtl and Mach numbers of gas. - The Prandtl number of the 
gas to use in equation (12) when determining the formula for vari- 
ation of the recovery factor Ag from experiments is based on the 

average static temperature T g ^ av . The Mach number M g ^ 3 , which 
has in the past generally led to good correlation of A g for static 

cascades, is calculated from the values of average relative velocity 
at the turbine outlet W 3 ^ aT , determined as described in the previ- 
ous section, and the static gas temperature T g ^ 3 . The formula is 

Mg,3 = - " 5 > a1r - (21) 

/ y 7 gS R sr r s J 3 

If surveys at station 3 are not feasible, the assumptions of 
' equal to rotor-blade-outlet angle and total conditions at sta- 

tion 3 equal to total conditions at station 4 provide means for 

• calculating V 3 in order to get W 3 and T g ^ 3 , which are needed 
to get T M g ^ av and M g ^ 3 . If the mixing losses are large, however, 
this procedure might lead to possible errors. 

Method of presenting Hq data. - For the case of experiments 
with cold air flowing around a static cascade of heated impulse 
blades, it is determined in reference 4 that the outside heat- 
transfer coefficients could be correlated for design angle of 
attack in the form of the following equation: 




16 


NACA BM E50A05 


where 2 0 /jc was the outside perimeter of the blade divided by it. 
The properties of the fluid, c p,g> M-gj and k g were determined 

at a film temperature that is equal to one-half the sum of the 
average blade -wall temperature and the static temperature of the 
fluid immediately upstream of the blade leading edges. The density 
and the velocity of the fluid, p g and V respectively, were 

determined from temperature and pressure measurements immediately 
upstream of the blade leading edge. The differences between the 
wall and stream temperatures were not very large. 

Further unpublished investigations have been conducted ori this 
cascade at design angle of attack using hot air (up to 300° F) and 
blade cooling. The difference between the wall and stream tem- 
peratures again was small. The data were worked up in the form 
given by equation (22) using film temperature, and so forth. The 
results when plotted gave a different line than the one obtained 
with cold air. Theory and experiments have indicated that if film 
temperature is used for correlating heat-transfer results, another 
factor in addition to Reynolds number and Prandtl number must be 
included in the equation. This factor is the ratio of the wall 
temperature to the static temperature of the fluid outside the fluid 
boundary layer on the surface (references 7 to 9). 

It is evident from experiments (references 8 and 9) and theory 
(reference 7) that if the properties of the fluid are based on the 
wall temperature rather than the film temperature and, in addition, 
according to reference 9, if the density p is also based on the 
wall temperature, correlation of data can be obtained over wide 
ranges of the ratio of wall to stream temperatures with the type 
of formula expressed by equation (22). This method was used for 
the case of hot- and cold-air experiments on the cascade of blades 
previously mentioned and a much better correlation of data was 
obtained than for the case where fluid properties were based on 
film temperature and the density was that determined from fluid 
measurements immediately upstream of the leading edges of the blades. 

A correlation of six sets of data obtained on cascades of 
turbine blades whose profiles were vastly different and with the 
flow at the design angle of attack has been attempted at the Lewis 
laboratory. The data of reference 10 indicated that using the wall 
temperature for properties including the density and a formula of 
the form of equation (22) was insufficient to obtain satisfactory 
correlation. It was found that a better correlation could be 
obtained if a formula of the type 
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Hu = f III 
(Pr) 1 ? 3 


Ke, 


(23) 


were used. A temperature- ratio effect is indicated "by theory 
(reference 11). It is therefore suggested that a type of formula 
like equation (23) he used in attempting to correlate the outside 
heat-transfer- coefficient data obtained in turbine experiments, 
or. 



vhere the properties of the fluid and the density are based on 
wall temperature. 

The foregoing attempts have been for design angle of attack 
of the gas flow. Changes in angle of attack of the fluid changes 
H q appreciably (reference 12). Thus the function f 17 in equa- 
tion (24) changes with each angle of attack and formulas like 
equation (24) must be established for each angle. 


The use of blade-wall temperatures to determine factors in 
the heat-transfer-coefficient equation necessitates locating 
thermocouples in the walls of the blades of the turbine. Usually 
such instrumentation is kept to a minimum but enough thermocouples 
should be used to give a fair average for the wall temperature. 


It is unfortunate that when analyzing the performance of the 
turbine for conditions other than those specifically tested, the 
blade-wall temperature must be assumed in order to calculate the 
heat-transfer coefficients from experimentally determined formulas. 
(The foregoing statement, of course, assumes that only by using 
blade temperatures can the data be correlated, which may not prove 
to be true in all future cases.) When the heat- transfer coefficient 
is calculated, the blade temperature distribution Is calculated 
according -to equation (3). If the calculated average blade tem- 
perature does not agree with the assumed value, then the procedure 
must be repeated using the new value. The results do not have to 
check exactly inasmuch as T-g^ av can vary somewhat without appre- 
ciably changing Eq. The heat- transfer coefficient also can vary 
somewhat without greatly affecting the results of equation (4). 
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In summary, for a cooled turbine it is suggested that experi- 
ments be conducted to establish a formula for each angle of attack, 
the form of which should be similar to equation (24), which it is 
expected will correlate the values of blade outside convection 
heat- transfer coefficients calculated according to equation (4)* 

The properties of the fluid, c -p f g> Hg> Pg> 831(1 k g should be 

based on the average blade temperature. The velocity Y g in the 

case of the rotor heat- transfer coefficients Is replaced with the 
average relative velocity. This average relative velocity Is equal 
to one-half the sum of Wg^ av and tiLe me ' fc k 0<is for deter- 

mining both of which have been given. The density should be based 
also on the average static pressure of the gas, which is one-half 
the sum of Pg,2 831(1 ^g,3* 

The static pressures of the combustion gas at the various 
stations is determined by instruments located radially and circum- 
ferentially so that a mass average can be found. 

The use of ijv. for the characteristic dimension in equa- 
tion (24) is suggested for the present. For turbines with low 
solidities, the blade may act like an Isolated airfoil; whereas 
with high solidities, a dimension such as hydraulic diameter of the 
passage between the blades will result In correlation of data from 
one turbine with those from another. In any one turbine, however, 
it is immaterial which dimension is used inasmuch as it is a con- 
stant and can be Included in the function f of equation (24) 
rather than In the Hu and Re terms. 


Inside Convection Heat-Transfer Coefficient 

Method of calculation. - The convection coefficient for the 
transfer of heat from the blades to the cooling air is calculated 
from data obtained In experiments on the heat gained by the cooling 
air, the average blade temperature, and the average effective air 
temperature in addition to the blade dimensions. 

For hollow blades such as those shown In figures 5(a) and 
5(b), the coefficient Is calculated using the formula 

Qe 

^ l i^ T B,av " T a,e,av) 


(25) 
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If an insert is used, such, as In figure 5(c), It is assumed that 
radiation is negligible and that the insert adds little heat to the 
cooling air. Consequently the temperature of the insert is neg- 
lected and equation (25) Is applicable to such a blade configuration. 

When fins are used in a hollow blade, as shown in figure 5(d), 
the heat-transfer coefficient E^ in equation (25) should be based 
on the average of the blade-vall and fin temperatures and the total 
surface In contact vith the air rather than on the average wall 
temperature and wall surface area, as in equation (25). Inasmuch 
as it is difficult to instrument turbine blades with thermocouples, 
it would be advantageous if thermocouples could be used only in the 
blade walls and that Hj_, based on average wall and fin temperatul-e 

and the total surface area, could be obtained from a coefficient 
ba3ed on wall conditions only. For. a finned blade such a relation 
exists. From reference 13, then, it can be 
blade the following relation is true: 

H* [~2 tanh (^Lg.) 

^ “ m + t L W~ 


shown that for the finned 


+ mj 


(26) 


where 



In equation (26), the coefficient Eg is based on wall conditions 
end is determined in a manner similar to equation (25 ) . 

It is necessary for correlation purposes to determine the 
coefficient H-j_ based on average wall and fin temperatures. Because 

it is difficult to solve for H-$_ in equation (26), it is advantageous 
to plot H-j_ against Ef and then for any value of Ef, determined 
from an equation similar to equation (25), can be found directly. 

Average effective cooling-air temperature T a,e,av* " The 

method of determining the average effective cooling-air temperature 
is very similar to that previously described for the effective gas 
temperature. The effective cooling-air temperature is defined as 
the adiabatic wall temperature of the blade and is related to the 
average total and static air temperatures by the expression 
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ipM 

“T 


a,av 

Ti — 

a,av 


a,e,av 

T 

x a,av 


(28) 


where K is the recovery coefficient of the "blade on the cooling- 
air side. 


As in the case of the recovery coefficient on the outside of 
the blade, a relation between cooling-air Mach and Prandtl numbers 
and Aa must be established from experiments conducted with con- 
ditions as nearly adiabatic as possible; that is, no gas is passed 
around the outside of the blades. The average blade temperature is 
inserted in equation (28) for to determine A a for these 

experiments. A relation of the form 


-p= - f V (Ma) (29) 

A ^ r a 


will probably correlate the data. The effective cooling-air tem- 
perature ^a,e,av> used to determine % from the heat- transfer 

experiments, is then calculated knowing the cooling-air conditions; 
namely, Pr a , Mg, T n a ^ a vJ T a ^ aV for these heat-transfer 

experiments and the relations given by equations (28) and (29). 
Consequently, in order to solve for T a ^ av , methods must be 

known for evaluating these four cooling-air conditions. 


Average static cooling-air temperature ' T a t aT . - It is desir- 
able to have temperature and pressure measurements at a number of 
spanwise stations along the blade in the cooling-air passages in 
order to have an accurate picture of the conditions of the cooling 
air. It Is assumed herein, however, that measurements are made 
only at the blade root and tip, which is usually more practical. 
Instruments should be so located In the blade passages at these 
two stations that accurate average measurements over the cross 
section can be obtained. Because the instruments are rotating with 
the blades, all total values calculated from the measurements are 
relative to the blades as previously mentioned. 


The effect of heat transfer, friction, and rotation on the 
’static- temperature distribution of the cooling air in the blade 
passage is such that the variation of this static temperature can 
be assumed linear if the Mach number of the cooling air is low, 
say below 0.5. Thus, 


T, 


T, 


a.av 


a.h + T < 


a,T 


2 


(30) 
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The static cooling-air temperatures at the root and. the tip, T a ^ ^ 
and T a ip, are determined by normal methods from the thermocouple 

measurements at these stations and the recovery factors of the 
particular instruments that are used. 

If the Mach number reaches high values, the static -temperature 
variation with the blade-length parameter y will probably be 
similar to that indicated in figure 6 and the assumption of linear 
variation is not quite valid. In order to accurately obtain the 
average static tempera ture of the cooling air under such conditions, 
a differential equation must be solved that involves the Mach numbers 
and the total temperatures at the root and the tip. The method of 
determining the values of total temperature mentioned have been 
described in connection with equation (10). The differential equa- 
tion and its solution to determine the variation of Mach number and 
total temperature through the cooling-air passages is the main sub- 
ject of reference 3. With the solution of the differential equation, 
the static temperature at any point x along the blade span is 
calculated from the expression 


rp« 


a,x 


a,x 


1 + 


^a - 1 


M a,: 


(31) 


The. average temperature T a ^ av is determined by integrating the 
values of T a ^ x with respect to the cooling-passage length, or 
blade span. 


Average total cooling-air temperature T" a av . - The total 

cooling-air temperature is practically linear vith respect to blade 
length for all Mach numbers and the average can be determined as 
the arithmetic mean of and T u a ^rp, which are the same 

values used in equation (10). 

Cooling-air Mach numbers , - The cooling-air Mach numbers at 
the blade root and tip can be calculated from the total and static 
pressures at these stations obtained from measurements by methods 
previously given. The correct Mach number to use in equation (29) 
to correlate the recovery-factor data is unknown at present but It 
is suggested that an average value for the cooling-air passages be 
used. For low Mach numbers, an arithmetic meaii of M a ^ and 

Mg^rp is thought to be sufficient. For high Mach numbers, an 

integrated average of the values determined by means of equation (l) 
of reference 3 is suggested. 
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Cooling- air Prandtl number. - The cooling-air Prandtl number 
to use in equation (29) should be “based on the average static 
cooling-air temperature T a ^ av . 

Method of presenting Hj_ data. - No data on the convection 
heat- transfer coefficient from “blades to cooling air in turbine- 
blade cooling passages using either cascades of blades or turbines 
are available. Consequently, no generalized form of a formula 
similar to equation (23) has been experimentally verified for these 
coefficients. The flov of cooling air through the turbine passages, 
however, can be thought of as flow of a fluid through a tube where 
the tube may have shapes such as those shown in figure 5. As a 
consequence, it would be expected that the inside coefficient data 
could be correlated by means of a formula similar in form to that 
which is known to correlate the experimental heat-transfer data for 
fluids flowing through tubes or pipes. The form of this formula for 
pipes is (reference 14, p. 164) 


Nu - f VI (Be,Pr) (32) 

A discussion is given in reference 1 on factors that may 
influence the heat-transfer coefficient of turbine-blade cooling 
passages. Effects of fluid flow at the passage inlet, which may 
cause to vary along the blade span, effect of an annulus shape, 

such as in figure 5(c), on the function f V3 - in equation (32), and 
other details are included. On the basis of this discussion, it is 
suggested that for the present, in accordance with equation (32), 
experiments should be conducted to determine if an equation similar 
to the following can be used to correlate the inside heat- transfer 
coefficients: 



The function f VI1 is determined from the investigations. 

It is suggested on the basis of pipe experiments that an expo- 
nent n of 0.4 be used in equation (33). The hydraulic diameter 

of the cooling passage to be used in equation (33) is defined 

by the table in figure 5 for the various blade-passage configurations. 
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It is suggested, on the "basis of results in reference 9, that the 
fluid properties, including density, in equation (33) "be based on 
average wall temperatures a v* Because, as discussed in refer- 
ence 1, analysis of data at inlet sections to pipes shows marked 
variation of heat-transfer coefficient with distance from the 
inlet, use of surface temperature for determining properties may 
not correlate the data for the coefficients in the blade passages 
where inlet effects may predominate. The use of surface temperatures 
was satisfactory for the long tubes used in reference 9 and is sug- 
gested as a first trial in attempts to correlate blade inside 
coefficients. 

The density should also be based on the average static pres- 
sure In the coolant passage. For low Mach numbers, this average 
can be the arithmetic mean of the pressures at the blade root and 
tip. The methods for determining the blade-tip pressures have been 
given. For high Mach numbers, the average static pressure must be 
obtained by Integrating calculated values determined as follows: 
Values of Mach number and static temperature are known through the 
passage by methods previously described. The relative velocity at 
each radial station is determined from 

w a,x ” *^a,x * l j7a.& R a T a,x ( 34r ) 

The static pressure at each radial station is then calculated 
from the expression 


•a.x 


w a B a T a,x 

w a>x A a,x 


(35) 


For low Mach numbers, the velocity W a to be used in equation (33) 
is an arithmetic mean of velocities at the blade root and tip deter- 
mined from measurements; for higher Mach numbers, W a is an inte- 
grated average of the velocities determined along the passage by 
means of equation (34) . All terms are thus determined so that the 
three parameters, Hu a , Pr a , and Re a of equation (33) can be 

calculated from the experimental data. 

The suggested form of formula as given by equation (33) Is 
based on the premise that the heat transfer from the blade to the 
coolant is due to forced convection. For heat transfer by forced 
convection alone, Reynolds and Prandtl numbers are the only sig- 
nificant parameters affecting heat transfer. Forced convection is 
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expected to be the predominant type of heat transfer for cooling air 
flowing through "blade passages that are open at both ends and are of 
the type used in air-cooled turbines. As a result, the use of equa- 
tion (33) is suggested as a first trial in attempts to correlate 
blade inside heat- transfer coefficients. 

Natural convection heat transfer has been studied through the 
use of water-cooled turbines. The coolant passages in the blades 
are holes that are plugged at the tip end of the blade. A circu- 
latory action of the coolant takes place, the coolant flowing from 
root to tip through the center of a hole and from tip to root along 
the outside of the same hole. A complete discussion of the govern- 
ing parameters for heat transfer for such a flow is given in refer- 
ence 1. With this type of circulation, the heat transfer is a 
natural- convection phenomena and in any such phenomena the Nusselt 
number is a function of Grashof and Prandtl numbers (references 14 
and 15). 

The possibility exists that, even in through coolant passages 
as in air-cooled blades, a small circulatory motion of the air in 
the passage may be superimposed on the main fluid motion, which is 
radial. This possibility may be more true for the case of hollow 
blades with large passages than for finned blades in which the 
passages are small tubes. Consequently, even with air-cooled blades 
of the form discussed (that is, with through passages), a small 
transfer of heat by natural convection in addition to a large trans- 
fer of heat by forced convection may occur. If the data do not 
correlate by a formula of the form of equation (33), then a form as 
follows may be applicable: 

Nu - f VIII (Gr,Ee,Pr) (36) 

Experiments must then be conducted over a range of both Grashof and 
Beynolds numbers to determine the function fVTII, 

The Grashof number for the turbine blade is calculated using 
the expression 



The gravity term g' is many times greater them, the ordinary 
gravity term g (reference 1) and varies with speed. This gravity 

term equals (0°v, where o> is the angular velocity of the 
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turbine -wheel and r is the radius from the center of rotation to 
the point considered. As a consequence g f varies from root to 
tip. It is suggested that g* be determined at a mean radius that 


is a geometric mean of the root and tip radii. 


<\/ r h r T* 


The coefficient of thermal expansion 0 a generally is 
expressed by the formula 


p _ dv/dT 
"^initial 


(38) 


where dv/dT is the change In specific volume with temperature and 
^initial is initial specific volume. Then equation (38) 
becomes for the cooling air In the blade passage. 


^a ~ Pa,h 

This equation reduces to 



a.h 


‘-■^a,h “ ^a,T 


s) 


T a,t - ^ Ta,T 


^a t 


'■i? 


a,h 


^■a,h ” -^a,T 


(39) 


(40) 


because, for small pressure changes, 3 is approximately equal to 
l/T. For small pressure changes in the turblqe cooling passage, 
that is, Pa,h/Pa,T s equation (40) reduces to £ a Z 

The methods of determining the air pressures and temperatures in 
equation (40) have been given. 


The density p a , the viscosity \x &t and the hydraulic diam- 
eter Dh,a in equation (37) can be treated In the same manner as 

suggested for the Eeynolds number determination (equation (33)). 


Miscellaneous Factors 

In giving the suggested methods for determining and correlating 
the outside and inside convection he at -transfer coefficients, Hq 
and from turbine investigations, methods for determining various 
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factors required to establish the validity of equation (3), which 
is part of any turbine- investigation program, have also been given. 
Methods have not been given, however, for obtaining Tg^ x and 

T a,e,h in equation (3) or various specific-heat and ratio-of- 
specific-heats terms in several equations. These terms are now 
given. 

Effective cooling-air temperature at blade root T ajQj jj. - 

The method of determining the effective cooling-air temperature at 
the blade root is the same as that for determining the average 
effective cooling-air temperature e ,av> except all values of 

Mach number, Prandtl number (used to get recovery factor at blade 
root), blade temperature, and so forth are the values obtained 
from the data at the blade root. The investigations for A a>ll 

can be simultaneously conducted with those for the average recovery 
factor inside the blade. The effective air temperature used in the 
equation to determine the variation of Aq ^ with Prandtl and Mach 
numbers is the integrated average of the peripheral blade tempera- 
tures at the root. 


Average blade temperature at any spanwise station 
The average blade temperature at any position x Tg^ x 


x 


l T B,x* “ 
is the 


integrated average of the blade-temperature readings at this posi- 
tion as obtained from the data. These values are compared with 
Tp x calculated using equation (3) to verify the use of this 
equation. Such averages have been mentioned as the first step 
required in the determining of T B,aV 


Average specific heat of cooling air in rotor blades. - An 
average specific heat of the cooling air is required in equation (10) 
to determine the heat gained by the air in the rotor-blade cooling 
passages. It is suggested that this specific heat be determined at 
an arithmetic average of the root and tip cooling-air static tem- 
peratures. This average should be accurate enough for all practical 
purposes . 


Average specific heat of cooling air in nozzle blades. - If 
the nozzle blades are cooled, the heat gained by the air passing 
through them is required. (See equation (15).) The specific heat 
in this equation can be determined in the same manner as suggested 
for for the rotor blades, except that temperature meas- 

urements made In the nozzle-blade cooling-air passages are used. 
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Over-all convection heat- transfer coefficient U. - A useful 
factor heretofore unmentioned is the over-all convection heat- 
transfer coefficient U. The heat transfer from the gas to the 
cooling air can he expressed in terms of the over-all coefficient 
as follows: 


Qa = m o bAT e (41) 

where AT e is the mean effective temperature difference between 
the gas and the cooling air. The over-all coefficient is useful 
in that if its variation with gas and coolant conditions is estab- 
lished, calculations of turbine performance at conditions other 
than those of the experiments are simplified. Also a check of the 
individual coefficients, Eq and H^, is obtained. Both of these 

objectives are obtained if the over-all coefficients calculated by- 
means of equation (41) agree with coefficients calculated from 
the formula 


_ 1 _ 

U = (42) 

HqIo + %*i 

where Hq and Hj_ are calculated from formulas established as 

previously described from experiments on the turbine (equations (4) 
and (25)). Equation (42) is valid with the assumption that radiation 
and temperature drop through the blade wall are negligible. It is 
recommended that values of U for all experiments made on a turbine 
be calculated by the two procedures and compared. Agreement will 
indicate that the assumptions of equation (42) are valid. If such 
agreement exists, it is also recommended that comparison of values 
of U using equation (42) with Hq and values calculated 

from the formulas that correlate the data, for example, equa- 
tion (24), be made with values calculated using the experimental 
data in equation (41) . Such a comparison will show the error 
involved in using Hq and Hi correlation formulas to calculate U 
for conditions other than those tested. If the error is small, U 
calculated in this manner can be used to calculate Qa for such 
conditions using equation (41) . The heat gained by the cooling air 
Qa Is required in performance evaluations. 

It may be possible to obtain a formula that correlates the 
over-all coefficients determined from the experimental data. 
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Because of the form of equation (42) and because Eq and Hi in 
turn are functions of several air and gas parameters, this formula 
may be complicated. If the values of U can be correlated, the 
steps required to determine for performance evaluations for 

other than experimental conditions vill be decreased. 

Mean temperature difference from gas to cooling-air AT e . - 

The over-all coefficient U is based on a temperature difference 
AT e (equation (41)), vhich is a mean of the temperature differences 
between the gas and the cooling air. This temperature difference is 
called the logarithmic mean temperature difference. The general sig- 
nificance of this factor is discussed in reference 14. 

In the case of the turbine rotor on the basis of equations (4) 
and (25), AT e is 


AT e 


i _ T , -j?— + 

g,e,av ■ L a,e,av 1^% 


(43) 


But from equation (41), 


AT e 


Q« 

UZ Q b 


Equation (42) can be obtained by equating equations (43) and (41). 

Charts have been made (for example, reference 16) by means of 
which AT e can be obtained for given inlet and outlet hot and cold 
fluid temperatures for many types of heat exchanger. As the air- 
cooled turbine can be considered to be roughly a cross-flow heat 
exchanger, it might be thought that the charts would be applicable. 
In deriving the equations for the charts, hovever, no distinction 
is made between effective, total, and static temperatures. The 
differential equations involve all such temperatures but each Is 
replaced by the fluid temperature. For the rotating turbine, such 
an assumption can cause large errors and It is recommended that 
AT e be determined from the difference between average effective 
gas and cooling-air temperatures, as shown by equation (43). 

Methods for determining these averages have been given. 
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S.tator-Blade Characteristics 

The foregoing discussion has dealt mostly with methods for 
determining the rotor-blade heat- transfer characteristics. If the 
nozzles or stator blades are to be cooled, experiments should be 
conducted to determine like heat-transfer characteristics for these 
parts . The methods are similar, the only difference being that 
velocities and so forth relative to the stator-blade surface are 
equal to the values relative to the turbine casing; consequently, 
in formulas suggested for the rotor blades, such symbols as W, 

T", and p" are replaced by Y, T* , and p* and the formulas 
can then be applied to the stator blades. The measurements outside 
and inside the stator blades are used to establish the functions of 
the formulas where required. 

Natural convection heat transfer (a Grashof number effect) is 
not considered in stator-blade characteristics. 


SUMMARY OF METHODS 

As a brief resume of the foregoing methods, the suggested pri- 
mary formulas for the determination of the rotor-blade heat-transfer 
characteristics of an air-cooled turbine are presented. The formulas 
with few changes, which are explained in the preceding section, 
apply equally well to the stator blades. 

In order to determine the coolant-flow requirements for safe 
operation at a given gas temperature for other than experimental 
conditions, the allowable and actual blade-temperature distribution 
must be calculated. The allowable- temperature curve is determined 
from stress-rupture data of the blade material, where the stress on 
a radial increment of the blade is calculated using equation (2). 

The actual blade-metal temperature distribution is found from equa- 
tion (3) if this equation is proved valid in the turbine experiments. 
In this equation, the convection heat-transfer coefficients on both 
the outside and inside surfaces of the blades must be known. For a 
given turbine, these coefficients are determined from the experi- 
mental data. 

The outside convection heat- transfer coefficient Is determined 
from the data using equation (4). After values of Eq have been 
evaluated, an attempt should be made to correlate these values In a 
dimensionless form similar to equation (24). It Is suggested that 
the fluid properties Including density be based on average blade- 
wall temperature. In addition, the density is based on the average 
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static pressure of the gas. The velocity W g in the Beynolds number 
is the average stream velocity in the "blade passages, relative to 
the blade. 

The inside coefficient Hj[ is treated in a similar manner and 
is determined from the experimental measurements using equation (25). 
Correlation of the values of so determined should be attempted 

by the use of dimensionless parameters in the form of equation (33) 
or equation (36), depending on vhich is the more applicable for the 
turbine investigated. As a first trial, the fluid properties, 
including density, in the parameters, are based on average blade 
temperatures. The density is also based on average static pressure 
in the coolant passages. The velocity V a is the average stream 

velocity relative to the blade and can be determined from an inte- 
grated average of values along the passage calculated by means of 
equation (34). The Grashof number in equation (36) is determined 
using equation (37). The coeff icient-of- thermal-expansion term in 
equation (37) is found by using equation (40). 

The purpose of evaluating and attempting to correlate the over- 
all heat-transfer coefficient U is to have a check on the individ- 
ual coefficients Eq and %, and also to simplify turbine- and 
engine-performance evaluations. The over-all coefficients are deter- 
mined from equations (41) and (42) and then compared to determine 
the error between Hq and Hi determinations and U determinations. 

As previously mentioned, it is necessary to know the effective 
fluid temperatures, both of the combustion gases and the cooling air. 
The recovery factor h, which is an indication of the effective 
temperature, is found from experiments by using the following 
formula: 


1-A = 


T H -T e 

mtt _m 


and can probably be correlated in a manner represented by 


- t (M) 

The mean temperature difference from the gas to cooling air 
AT e was used in the determination of the over-all heat-transfer 

coefficient in equation (41) and is determined from equation (43). 


A 

V pF 
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The heat loss to the cooling air is the term on which the 
convection heat- transfer coefficients are directly dependent and 
can he calculated hy using equation (10). The factor f s is deter- 
mined by motoring the turbine wheel at several speeds with no gas 
passing across the blades and with various quantities of cooling 
air passing through the blades. In an experiment of this type 
qg^ = 0, and f s can be readily determined. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio . 
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APPENDIX A 
SYMBOLS 


A 

B 

b 

C 

c p 

D 

D h 


The following symbols are used in this report: 
cross-sectional area, sq ft 
force, lb 

height of blade (span), ft 

nozzle-discharge coefficient 

specific heat at constant pressure, 
Btu/(lb)(°F) 

characteristic dimension in Eeynolds and 
Nusselt numbers, ft 

hydraulic diameter 4A/Z , ft 

slip factor 


fl, f II, fill, f I7, 

unknown functions 

t v , f VI , fVII.fVIIX 


Gr 


Grashof number 


D 5 p 2 g«p( T B -T e ) 


g 


g* 


ratio of absolute to gravitational unit of 
mass, Ib/slug, or acceleration due to 

o 

gravity, ft/sec c 

2 

acceleration used in Grashof number ((0 r), 
ft/sec^ 


H 

J 


convection heat-transfer coefficient, 
Btu/(sec)(sq ft)(°F) 

mechanical equivalent of heat, 778.3 ft-lb/Btu 


k 


thermal conductivity, Btu/(°F)(ft)(sec) 
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L f 

l 

M 

m 

Nu 

n 

Pr 

P 

P* 

p" 

Q 

q 

E 


effective width, of fins from hlade vail to 
mean camber line, ft 

perimeter, ft 

Mach number 

average spacing between fins in finned 
blade, ft 

Nusselt number, 
exponent 

Prandtl number, (¥) 

static pressure, lb/sq ft absolute 

total pressure relative to turbine casing, 
lb/sq ft absolute 

total pressure relative to moving blade, 
lb/sq ft absolute 

heat flow, Btu/sec 

heat flow per unit weight, Btu/lb 

gas constant, ft-lb/(lb) (°F) 



Ee 

r 

S 

T 

^ e 
T* 


Beynolds number, 

radius, ft 
stress, lb/sq in. 
static temperature, °B 

mean temperature difference effecting heat 
transfer, °F 

total temperature relative to turbine casing, °E 
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total temperature relative to moving blade, °E 

over-all beat- transfer coefficient, 
Btu/(sec)(sq ft) (°F ) 

tangential velocity, ft/eec 

absolute velocity, ft/sec 

specific volume, cu ft/lb 

relative velocity, ft/sec 

weight flow, lb/sec 

W 2 r ll w a 

J S H o^o^g,e “ T a,e,h) 

radial distance from blade root to point on 
blade span considered, ft 

(0 2 v a 2 (X + 1) c Pt a 

JgHo 1 0 (^g,e - ^a,e,h} 


b 


1 

XTT 


HqZ o b 

w c 
a p,a 


x 

b 


angle between absolute and tangential velocity 
vector, deg 

coefficient of thermal expansion, l/°F 


ratio of specific heats. 


( °p/ °v) 


recovery factor 
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P 

t' 

0 

n 

CD 

Subscripts: 

a 

at 

av 

B 

e 

f 

e 

h 

i 

o 

S 

T 

t 

x 


absolute viscosity, slugs/(sec) (ft) 
density, slugs/cu ft 

thickness of fins in finned blades, ft 



V k B T 

vork added to cooling air by blade rotation, 
Btu/(lb) 

angular velocity, radians/sec 

cooling air 
atmospheric 
average 
blade 

effective, used vith symbol for temperature 
and denotes temperature effecting « heat 
transfer 

finned blades 

combustion gas 

blade root 

inside surface of blade 

outside surface of blade 

stator blades 

blade tip 

pressure tube 

point along blade span 
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1 station in combustion- gas stream at turbine 

inlet 

2 station in combustion- gas stream between 

nozzles and rotor blades 

3 station in combustion gas stream at rotor 

outlet 

station in combustion gas stream downstream 
of rotor where swirl has largely disappeared 


4 
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APPENDIX B 

CORRECTION OF RECORDED MANOMETER PRESSURES OBTAINED 
WITH ROTATING TUBES 

Consider the pressure equilibrium of a fluid element in a 
pressure line from the center line of the turbine to some radial 
point in the blade coolant passage (fig. 7). The pressure force 
directed outward on the fluid element is 

B 1 “ Pa-H ( B1 ) 


where 

p a pressure of cooling air in tube, lb/sq ft absolute 
A-t cross-sectional area of tube, sq ft 

The pressure force directed toward the center of rotation is 

B 2 * ( B a + ^ Ar t ) At (BE) 


where 

dp a 

— - differential change of pressure in tube with respect to 

t change of radial length, (lb/sq ft)/ft 

Ar^ radial increment of fluid considered, ft 

The difference between these two forces gives a resultant 
inward force of ( dpg/tir^. ) Ar.|-A^ , which must be balanced by the 

centrifugal force on the element in order for equilibrium to 
exist. This centrifugal force Bj can be given by the formula 


b 3 81 [fa,t A t &r t] ® 2 r t 


where 

P a ,t density of fluid in tube, slugs /cu ft 
« angular velocity of tube, radians/sec 


(B3) 



38 


NACA RM E50A05 


r t radius to midpoint of element considered, ft 

Because B 2 equals B 1 plus B 3 , the following equation results 
from equations (Bl) to (B3): 


d P a 2 

dr t = p a,t w r t 


(B4) 


When the fluid is compressible and the angular velocity is large, 
the density must be considered variable with pressure and tempera- 
ture. Combining equation (B4) with the equation of state for a 
perfect gas gives 


_ — Pa r 

SE a T a ,t * 


(B5) 


where 

R a gas constant for air in tube, ft-lb/(lb) (°F) 

* 

o 

temperature of air in tube, R 

The temperature of the air will probably not be constant 

along the tube. The selection of the proper temperature is important 
inasmuch as a 100° F difference between the correct average tempera- 
ture and the temperature used can in some cases cause an error in 
pressure of more than 1 pound per square inch. Thermocouples should 
be inserted in the pressure tubes over the length of the tube to 
find this average temperature. It is possible that the tubes if 
fastened to the rotor disk will assume the temperature of the disk 
and an average temperature of the disk will be adequate. 

Equation (B5) can be put in the form 

dPa 

Pa 

which after integration becomes 

2 2 

r t . , 

a * 2 6Va,t 


0) 


S^a T a,t 


■ t art 


l°6e P, 


(B6) 
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The integration constant K must "be evaluated. When r t = 0, 

K equals the logarithm of the pressure at the center line of the 
turbine, which is the value measured on the manometer board. This 
pressure is Pa,m> where m indicates measured. Then equation (B6) 
becomes 


l°g e 


2 2 
0) r t 


Pa,m ^S^a T a,t 


or 


Pa 

Pa,m 


e 


' 2 2 \ 
m r t ) 
2gR a T a ^t/ 


The manometer-board reading is the difference between the pressure 
p a and the pressure of the room in which the manometer board is 

placed or 

A Pa,m = Pa,m ^ Pat 

The true pressure at any radius r^ in the blade passage is then 


Pa " ( A Pa,m ± Pat) e 


2sR a T a,t 
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air flow 


Centrifugal force + B 



Cutaway view of typical 
turbine blade. 


Area Ag at base of Increment 

Cross-sectional area Ag av at 
center of Increment 


(b| Increment of turbine blade used 
in cen t r l f uga I -s t ress calculations. 


Figure 2. - Typical air-cooled blade with Increment to be used In 
centrifugal-stress calculations. 



Pressure-ratio function, ( — ^ 


44 


NACA RM E50A05 


. oa 



Figure 3. - Relation of pressu re- rat i o function, 
pressure ratio, and temperature ratio for several 
values of ratio of specific heats of gas. 


Temperature ratio. 
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Figure 4. — Typical velocity diagram for a turbine blade. 
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Tubes 


Ho I f ow blade 


Blade with Insert 


F I nned blade 


a s s a g e area 


.wetted perimeter 


4 x passage area 
sum of Inside blade 
perimeter and Insert perimeter 


4 x summation of cross- 
sectional area of each passage. 


summation of Individual 
passage perimeters 


I n se r t 


Figure 5. - Four typical air-cooled turbine-blade configurations 
and hydraulic diameter of cooling-air passage to be used In 







CooIFng-alr temperature, 


NACA RM E50A05 


47 



Figure 6. — Variation of Mach number, total temperature, 
and static temperature of cooling air from blade 
tip to root . 
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